A repulsive atomic gas in a harmonic trap on the border of itinerant ferromagnetism 
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Alongside superfluidity, itinerant (Stoner) ferromagnetism remains one of the most well- 
characterized phases of correlated Fermi systems. A recent experiment has reported the first evi- 
dence for novel phase behavior on the repulsive side of the Feshbach resonance in a two-component 
ultracold Fermi gas. By adapting recent theoretical studies to the atomic trap geometry, we show 
that an adiabatic ferromagnetic transition would take place at a weaker interaction strength than is 
observed in experiment. This discrepancy motivates a simple non-equilibrium theory that takes ac- 
count of the dynamics of magnetic defects and three-body losses. The formalism developed displays 
good quantitative agreement with experiment. 
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The transition between ferromagnetism and param- 
agnetism with increasing temperature is considered to 
be a canonical example of a continuous phase transi- 
tion, and the phenomenon has been firmly established in 
many materials. The ability to tune the critical temper- 
ature through varying pressure presents experimentalists 
with an opportunity to investigate the novel many-body 
physics predicted to arise in the vicinity of a quantum 
critical point However, whether this phase behav- 
ior derives from soft magnetic fluctuations, which are be- 
lieved to drive the Stoner transition first order [H, 0, @ , or 
is a consequence of coupling to auxiliary degrees of free- 
dom such as lattice vibrations remains a subject of con- 
siderable interest and debate [f| . The control afforded by 
Feshbach resonance phenomena in ultracold degenerate 
Fermi gases presents a controlled platform from which to 
explore strongly correlated repulsive Fermi systems, in- 
cluding collective effects from non-equilibrium dynamics 
and spin textures [lol , [ll[ , and quantum critical phenom- 
ena 

Although a few early experiments @, [1] conducted on 
the repulsive side of the resonance hinted at ferromag- 
netic behavior Q , these investigations were hindered by 
the challenges posed by the cold atomic gas set-up, with 
the fixed relative populations of particles, trap confine- 
ment, atom loss through three-body interactions, and 
non-equilibrium physics, rendering the conclusive iden- 
tification of ferromagnetism impossible. However, in a 
recent study, Jo et al. Q succeeded in observing the first 
strong evidence for novel phase behavior, consistent with 
itinerant ferromagnetism in an atomic gas of ~ 6.5 x 10 5 
6 Li atoms. To overcome the obstacle of atom loss through 
three-body interactions the experiment was carried out 
under non-adiabatic conditions, with the atoms prepared 
in the disordered non-ferromagnetic state and the mag- 
netic field ramped to the repulsive side of the resonance 
in 4.5ms and then held fixed for a further 2ms. To assess 
the viability of the experimental design, it is crucial to 
have detailed predictions for the expected phase behavior 
in the trap geometry. In the following, we will adopt two 
stands of investigation. Firstly, it is important to under- 
stand what we expect to see in an idealized equilibrium 



trap geometry. This study provides a benchmark both to 
assess the current experiment and to guide future studies. 
However, while this analysis achieves a qualitative agree- 
ment with experiment, it also exposes important discrep- 
ancies highlighting the need to consider non-equilibrium 
effects. In the second part of the letter we turn to ad- 
dress non-equilibrium dynamics to critically analyze the 
experimental observations and establish strong quantita- 
tive agreement between theory and experiment. 

An atomic gas at equilibrium provides an ideal plat- 
form from which to analyze experiment. For a two- 
component Fermi gas with a local contact interaction 
and equal masses, the Stoner Hamiltonian is defined by 
H = J d 3 x [^2 a a^p 2 /2m)a a + ga]-a^aia^], where the 
atoms are distinguished by a pseudospin index a G {t> !}• 
Taking into account all contributions to second order in 
g, the energy density of the Stoner Hamiltonian for a 
spatially uniform system can be expressed as @, 0, H[ , 
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where n^(e) = 1/[1 + c' 3 '- 6- '^'] denotes the Fermi dis- 
tribution, with spectrum et, chemical potential fi a , and 
reduced temperature (3, N a is the number of particles, 
v is the density of states at the Fermi surface, fcp is the 
Fermi wave vector of the corresponding non-interacting 
system, and the scattering length a fully characterizes the 
strength of the contact interaction close to resonance. We 
also set h = m = 1. Retaining only the leading interac- 
tion correction, 0(kpa), Eqn.[T] recovers the conventional 
mean-field Stoner theory, which predicts that the tran- 
sition remains second order down to zero temperature. 
However, the soft transverse magnetic fluctuations [l2j 
encoded in the second order term [l3| have the effect of 
driving the ferromagnetic transition first order 0, 0, Q 
at low temperature. 

To address the non-uniform atomic trap geometry, we 
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must determine the free energy density, / = e — /z[n-|-(r) + 
n I (r)] . Here the Lagrange multiplier fj. enforces the con- 
straint of a fixed total particle number imposed by the 
trap geometry. It is also convenient to rotate the basis to 
the axis of net magnetization s 6 {+,—}. To analyze the 
spatially inhomogeneous atom distribution we invoke a 
local density approximation which allows the variational 



minimization 



U 
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This leads to two equations for the 



effective local chemical potentials //±(r) of the majority 
and minority species in the rotated basis that must be 
solved self-consistently, 
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within the external confining potential V(r) = wr 2 /2. 
To provide a reference with experiment 0], we describe 
the interactions in terms of a dimensionless interaction 
strength expressed in terms of the Fermi wave vector k F a 
at the center of a trapped non-interacting Fermi gas. 

We first address the behavior of the cloud size across 
the range of interaction strengths. Fig. |Tfb) confirms 
that, for weak interactions, the RMS cloud size follows 
the universal scaling relationship i?Q MS (l + g^" 2 ^ k F a) 
that can be derived from Eqn. [2] However, fluctuation 
corrections soon increase the cloud pressure causing it 
to inflate rapidly to its fully polarized value 2 1 / 6 i?Q MS . 
As cloud density is suppressed by pressure, the onset of 
ferromagnetism at the center of the trap takes place at 
an enhanced interaction strength compared with the uni- 
form system. Secondly, in Fig. [TJc) we study the total 
kinetic energy of the atoms, which is in agreement with 
the mean- field prediction of LeBlanc et al. With 
increasing interaction strength the local density falls so 
that the kinetic energy decreases to §-Ep(l — 945^ hp a). 
The fluctuation corrections increase the cloud pressure 
and reduce the density so driving the kinetic energy 
downwards still further. The release energy recorded 
by experiment extrapolated to zero interaction strength 
is ~ 0.45E"p, which is marginally higher than the an- 
alytical prediction of 0.375i? F . This discrepancy could 
arise through the recovery of interaction energy upon 
release from the trap, remnant eddy currents in the 
coils, or non-equilibrium effects. The final experimen- 
tal probe that we address here is the atom loss rate due 
to three-body recombination. Integrating the loss rate 
r = T (k F a) 6 J n+(r)n_(r)[n + (r) + n_(r)]d 3 r [U] over 
the trap yields the variation with interaction strength 
shown in Fig. [ljd). The mean-field case agrees with 
the prediction of LeBlanc et al. [llj. The initial rapid 
rise in the loss rate can be attributed to the increase 
in the (k F a) 6 coefficient, through the universal rela- 
tion r/r s» 7(fc F a) 6 /i 6 /2887r 4 w 3 / 2 . The later decay in 
the loss rate is due to the suppression of the product 
n+(r)n_(r) with the onset of magnetization. 
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FIG. 1: (Color online) (a) Atoms remaining after the hold 
time, (b) RMS cloud size, (c) release energy, and (d) loss rate 
with the dimensionless interaction parameter k F a shown on 
the primary x-axis in the mean-field case (green line), with 
fluctuation corrections (red line), and with defects (magenta 
line). The dashed blue line shows the trend in the small kpa 
limit, the cyan line is the Stoner criterion, and the experi- 
mental points of Ref. || are also highlighted. 



Although there is reasonable qualitative agreement be- 
tween theory and experiment, a marked divergence arises 
in the experimental prediction of the interaction strength 
at the onset of ferromagnetism at kpa « 2.2. The the- 
oretical prediction from mean-field theory is kpa w 1.9, 
whereas it is at kpa « 1.1 if fluctuation corrections are 
taken into account. This discrepancy prevents us from 
drawing a definitive conclusion on whether the transi- 
tion is first or second order. However, in future ex- 
perimental studies, it should be possible to gain insight 
into this question by studying the momentum distri- 
bution of the atoms shown in Fig. [51 The latter can 
be directly measured by studying the spatial distribu- 
tion following a ballistic expansion (3]. A single non- 
interacting atomic species has the momentum distribu- 
tion rik oc k 2 (2/i — fc 2 ) 3 / 2 . In the partially polarized 
regime the atomic gas will be ferromagnetic in the high 
density regions - at the center of the trap. Since the 
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FIG. 2: (Color online) The momentum distribution n^./N of 
a half net polarized atomic gas in (a) the mean-field case and 
(b) with fluctuation corrections. The cyan filled curve shows 
the population of majority spin particles and the blue filled 
curve shows the minority spin particles. The dotted curves 
show the expected distribution for the majority spin particles 
in a harmonic trap, oc k 2 (2(i — k 2 ) 3 ' 2 , and simply due 
to the density of states, rik oc k 2 . The inset curves show 
the radial distribution of atoms (orange, primary y-axis) and 
magnetization (blue, secondary y-axis). 



transition from an unpolarized to a fully polarized gas 
takes place over a small range of interaction strengths, 
the atomic gas will be partially polarized over a small 
range of radii. The remainder of the gas is unpolarized. 
The thin shell of partially polarized gas gives rise to a 
sharp double peak structure in the momentum distribu- 
tion in both the mean-field case and also when fluctuation 
corrections are taken into account. However, the mean- 
field ferromagnetic transition takes place over a larger 
range of radii than the fluctuation corrected case, which 
smears out the double peak feature. 

To assess whether the observed experimental phe- 
nomenology [i| is consistent with the development of 
itinerant ferromagnetism, it is important to understand 
the source of the discrepancy in the predicted interaction 
strength. In doing so, we will exploit the marginally adi- 
abatic nature of the experiment to explore the quench 
dynamics of the transition. This necessitates taking into 
account two contributing factors: how the ferromagnetic 
state condenses out of the normal phase, and the renor- 
malization of the interaction strength due to three-body 
losses. The emergence of the ferromagnetic state is non- 
trivial since the propagation of the condensed spin align- 
ment is bounded by the spin wave velocity. The quench 
propels the system deep into the ferromagnetic regime, 
so monopole defects are condensed out from the param- 
agnetic phase [H, Ell • As the defects grow they mutu- 
ally annihilate, delaying the formation of the ferromag- 



netic phase, which in turn increases the effective inter- 
action strength required to observe ferromagnetic phe- 
nomena. The second component of our dynamical study 
are three-body interactions. These reduce the atom den- 
sity and therefore abate the effective interaction strength 
fcpa, which to compensate must be artificially raised to 
observe ferromagnetism. The interplay between these 
two dynamical effects is expected to raise the interaction 
strength required to observe ferromagnetism to be more 
in accord with the experimental findings. We now detail 
how to separately incorporate these effects into a descrip- 
tion of the atomic gas, before merging them into a single 
formalism. Since the characteristic timescale associated 
with temperature, 0.5ms, is similar to time scales of the 
Ketterle experiment @, we can treat the atomic gas as 
if it were at zero temperature. Due to the short time- 
scale of the quench, for simplicity we can focus attention 
on just the mean-field component of the theory [l7j and 
make the approximation that the ferromagnetic transi- 
tion takes place at the Stoner criterion fcpacRiT = 7r/2. 

The ferromagnetic quench deep beyond the spinodal 
line leads to the condensation of topological monopole 
defects (T^. In each defect the spins radiate out from a 
central core with magnetization <f>(r > £) = 4>o(l—£ 2 /r 2 ), 
where £ = l/kFy/kpa/kpacmT — 1 is the healing length. 
A uniform atomic gas would be only partially polar- 
ized over a range of interaction strengths A w 0.29, 
so therefore during the field ramp over time i ram p to 
a target interaction strength /cfoaim, the gas will be 
partially polarized for the time t iamp A/kpaAiM- The 
study of Babadi [TtJ shows that after this characteristic 
quench time the defects will adopt an initial size Lq = 

[(fi,/77ifc^ 3 ^ 2 /67r 2 )t ra mpA/fcFaAiM] 1 ^ 3 - Subsequently, the 
defects grow and annihilate as they compete to select the 
equilibrium spin. This expansion of the defects can be 
summarized by the growth law of yjht/m (l9j . Combined 
with the initial condensed size of the defect s, we there- 
fore model the defect size by L(t) = X(L a + yjht/m), with 
A being an unknown dimcnsionlcss constant that we in- 
clude to account for the system dependent dynamics that 
go beyond the current analysis. 



Having captured defect growth in the model, we now 
incorporate it into a formalism that also accounts for the 
atom loss due to three-body interactions. We employ the 
standard atom loss rate formula [l5| to describe three 
body recombination h = — lll(na 3 ) 2 mjn^/7m, where 
e = 4.56/i 2 n 2 / 3 /m is the average kinetic energy. The fall 
in loss rate due to Pauli blocking is described by a geo- 
metric term that can be expressed in terms of the mag- 
netization 4> through n^ni/n 2 = (n 2 — cf> 2 )/An 2 . As the 
spins within a defect are not parallel, three-body recom- 
bination occurs within a defect, which can be expressed 
through the geometric term as n^n^/n 2 = £ 2 /2r 2 +0(£ 4 ). 
Integrating this loss over one defect predicts a net loss 
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rate of 

3 dt K 1 h J L(0) 4r2 

^ = -0.144-n 5 / 3 (fc F a) 6 [L(<)-L(0)]C 2 L- 3 (t).(3) 
dt m 

Starting in the disordered regime, following the magnetic 
field ramp we numerically propagate the atom loss rate 
forwards in time. If the system breaks into the ferromag- 
netic regime, L(t) starts to grow as the defects expand 
and annihilate. Should atom loss cause the system to 
reenter the normal regime we immediately revert to the 
L = limit. To determine the atomic distribution within 
the trap we self-consistently evaluate the local effective 
potential following the method employed for the uniform 
case in Eqn. [2] 

In Fig. [TJa) we first analyze the population of rem- 
nant atoms following the field ramp. At weak interaction 
strengths the number of atoms lost to three body recom- 
bination scales as ~ (fc F a) 6 . The onset of ferromagnetism 
reduces the atom loss, as the more rapid onset of ferro- 
magnetism and smaller core length conspire to cut losses. 
Deep within the ferromagnetic regime, around 80% of the 
atoms remained in the experiment Q , which was used to 
calibrate the constant A = 5. Having determined the di- 
mcnsionless constant, we now monitor the ramifications 
of this choice on the growth of the defects. The ratio 
of the defect size to the core radius, L/£, after the hold 
time, increases with fcpa due to the proliferation of de- 
fects commencing earlier and the falling core radius. The 
rise of L leads to a fall in the number of defects remaining 
after the field hold to ~ 8. However, no firm evidence has 
yet been reported for domain formation in experiment. 

Now that we have calibrated the theory against the 
experimental results we can now proceed to compare the 
other physical observables against the experimental data. 
We first look at the RMS cloud size in Fig. QJb) . In the 
weakly repulsive regime the cloud size increases with in- 
teraction strength due to the enhanced pressure between 
atoms. However, at interaction strengths 1 < k^a < 2.2 
the increasing atom loss dominates, causing the cloud to 



shrink. The characteristic minimum formed is in accord 
with the experimental results, and critically is delayed 
until the interaction strength fcjja ~ 2.2 both because the 
loss of atoms hinders the onset of ferromagnetism, and 
since on the border of ferromagnetism the longer core 
length increases atom loss. Upon ferromagnetic ordering 
the atom loss falls off causing the cloud to re-expand, 
which eventually exceeds its original radius. With weak 
interactions the kinetic energy in Fig. [TJc) falls both as 
atoms are lost and the cloud dilates. Upon reaching the 
critical interaction strength of fcpa f=a 2.2 ferromagnetic 
ordering takes place and the kinetic energy rises as the 
ordered atoms posses a larger Fermi surface. Across the 
whole range of interaction strengths the theoretical pre- 
diction for the kinetic energy falls short of the experimen- 
tal results, which could be due to the gas being inherently 
out of equilibrium, recovery of some interaction energy 
upon release from the trap, or remnant eddy currents in 
the coils. At high interaction strength the kinetic energy 
continues to rise as the atom loss decreases, indicative 
of the upturn seen in the experimental results. Finally, 
in Fig. [IJd) we address the atom loss rate which climbs 
strongly as ~ (fcpa) 6 before the onset of ferromagnetism. 
Consistent with the other experimental probes, the peak 
atom loss is delayed until fc F a w 2.2. Following the onset 
of ferromagnetism the atom loss drops away more slowly 
than without dynamic effects due to the necessity for the 
atom loss rate to grow and the initial atom loss during 
the ramp time. 

In conclusion we have performed a detailed critique of 
the results of the first experimental signs of ferromag- 
netism in a cold atom gas. Firstly, we demonstrated that 
the results are not consistent with the formation of an 
equilibrium ferromagnetic phase. Secondly, we demon- 
strated how the experiment indicates that the ferromag- 
netic phase is formed through the condensation of defect- 
antidefect pairs that subsequently undergo mutual anni- 
hilation, which delays the formation of the ferromagnetic 
phase to an enhanced interaction strength of k-pa = 2.2. 

We thank Eugene Demler, Jonathan Edge, Andrew 
Green, Zoran Hadzibabic, and especially Gyu-Boong Jo 
and Wolfgang Kettcrlc for useful discussions. 



* Electronic address: gjc29@cam.ac.uk 
[1] J. A. Hertz. Phys. Rev. B 14, 1165 (1976); A.J. Millis. 

Phys. Rev. B 48, 7183 (1993). 
[2] D. Belitz, T.R. Kirkpatrick and T. Vojta. Rev. Mod. 

Phys. 77, 579 (2005). 
[3] A. A. Abrikosov and I.M. Khalatnikov. Soviet Phys. 

JETP 6, 888 (1958); F. Mohling. Phys. Rev. 122, 1062 

(1961). 

[4] R.A. Duine and A.H. MacDonald. Phys. Rev. Lett. 95, 
230403 (2005). 

[5] G.J. Conduit and B.D. Simons. Phys. Rev. A 79, 053606 
(2009). 

[6] R.B. Laughlin et al, Adv. Phys. 50, 361 (2001). 
[7] S. Gupta et al, Science 300, 1723 (2003). 



[8] T. Bourdel et al, Phys. Rev. Lett 91, 020402 (2003). 
[9] G.-B. Jo et al, Science 325, 1521 (2009). 
[10] I. Berdnikov, P. Coleman and S.H. Simon. Phys. Rev. B 

79, 224403 (2009) 
[11] L.J. LeBlanc et al, Phys. Rev. A 80, 013607 (2009). 
[12] D. Belitz et al, Phys. Rev. B 58, 14155 (1998). M. 

Shimizu. Proc. Phys. Soc. 84, 397 (1964); D. Belitz, T.R. 

Kirkpatrick, and T. Vojta. Phys. Rev. B 55, 9452 (1997); 

J. Betouras, D. Efremov and A. Chubukov. Phys. Rev. B 

72, 115112 (2005); D.V. Efremov, J.J. Betouras and A. 

Chubukov. Phys. Rev. B 77, 220401(R) (2008). 
[13] G.J. Conduit, A.G Green and B.D. Simons. 

larXiv:0906.1347l 

[14] W. Ketterle and M.W. Zwierlein. larXiv:0801.2500l 



5 



[15] D.S. Petrov. Phys. Rev. A 67, 010703(R) (2003). 
[16] G.-B. Jo and W. Ketterle (private communication). 
[17] M. Babadi, D. Pekker, R. Sensarma, A. Georges and E. 
Dernier. larXiv:0908.3483l 



[18] A.J. Bray. Adv. Phys. 43, 357 (1994). 
[19] J.S. Langer, Solids far from equilibrium, edited by C. 
Godreche (Cambridge University Press) (1992). 



